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Abstract 

Stochastic neutron transport theory is applied to the derivation of the two-neutron- 
detectors cross power spectral density for subcritical assemblies when external pulsed 
sources are used. A general relationship between the two-detector probability gene- 
rating functions of the kernel and the source is obtained considering the contribution 
to detectors statistics of both the pulsed source and the intrinsic neutron source. An 
expansion in a-eigenvalues is derived for the final solution, which permits to take 
into account the effect of higher harmonics in subcritical systems. Further, expres- 
sions corresponding to the fundamental mode approximation are compared with 
recent results from experiments performed under the MUSE-4 European research 
project. 



1 Introduction 



In last years, researchers have shown an increasing interest on the conceptual 
development of accelerator-driven systems (ADS) for nuclear waste transmuta- 
tion and energy production purposes. An important issue regarding its future 
industrial applicability is the development of a periodically subcriticality level 
measurement and monitoring technique, since both its operation safety and 
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its performance as a part of the nuclear fuel cycle shall be seriously affected 
by this variable. 



Following the study of the neutron fluctuations in a multiplying medium, 
Courant and Wallace! ( 1947 ). several static and dynamic methods have been 
proposed and studied for years concerning their appli cability for th e determi- 
nation of some nuclea r reactor physics parameters f see lUhri 3 (|l97nh : lw7lliamsl 
f)l974f) : iLewinsI (|l978h : ICarta and D'Angelol (|l999h . Within the group of dy- 
namic techniques, apparently the utilisation of the neutron-fission chain fluctu- 
ations for nuclear assemblies subcriticality determination was firstly suggested 
by Bruno Rossi (the Rossi-a method). In these methods neutron detector 
counting rates related to individual fission-chain events must be discerned 
from the total counting rate, therefore these m ethods are applicable to sub - 
critical systems near delayed critical condition, I fWta acd H^wS B 



Further, dynamic methods based on a time-dependent external neutron source 
were proposed bv lPerez et al\ (J1964J). Recently, the use of these methods has 
increased during the MUSE European experimental studies carried out at the 
MASURCA facility (Cadarache, France) due to their applicability in order to 
investigate ADS kinetic parameters. In these experiments D-D and D-T neu- 
tron sources running in pulsed mode have been used, although the utilisation 
of an external spallation proton source has also been thought. An y way, mea- 
surem e nts can be done in t wo dif ferent ways ([Valentine et al. iDegwekerl 
f)2003h : ICeder and Pazsitl (fcOOSl )): in the first one, the neutron detector time 
gate is synchronized with the external neutron pulse injection, thus this me- 
thod is referred to as deterministic pulsing method, whereas, in the second 
case, the relative delay between the neutron pulse injection and the begin- 
ning of the neutron counting time is uniformly sampled between zero and the 
pulsed source period, which is known as stochastic pulsing method. For the 
latter case, the neutron source can be assumed to have the form 



S(t)=k £ 5(t-(S + mT)) 



where k is the number of protons injected per proton pulse, T is the pulsed 
source period, and £ is uniformly sampled within the time interval [0, T]. 
Obviously, for the deterministic pulsing method we will put £ = 0. 



On the other hand, IPall ( 19581 ) developed a general theory for the stu dy of 



the st ochastic neutron field. This model, complemented afterwards by Belli 



(1965). completely describe the stochastic neutron field in a fissile assembly. 



Later, iMufioz-Cobo et al\ (|1987I ) derived expressions for the variance of the 
number of counts in a detector and the CPSD for a Poissonian source without 
delayed neutrons from the general neutron stochastic transport theory. These 
approaches permit to go beyond classical point kinetic approximations, taking 
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into consideration general problems with spatial, spectral, and angular depen- 
dence. They have been extensively applied to nuclear subcriticality safety and 
non-destructive nuclear fuel assay problems. 

Classical reactor noise methods are no t correct when pulsed o r corr elated 
sources are used ( Matthes et al. ( 1988h : Behringer and Wvdler For 



these kind of sources, forward Kolmogorov's approach is not valid because 
of the non-Markovian character of the process, while back ward Green ' s func - 



tion description needs to go beyo nd Poissonian behaviour ([Degweker (2003) 
Ballester and Munoz-Cobol ifeOOSh ). 



In our work we have derived a relationship between the source probability 
generating function and the kernel probability generating function when non- 
Poissonian spallation neutron sources are considered. We have also studied 
the effect of the intrinsic neutron source due to spontaneous fission occur- 
ring in major actinides forming part of the nuclear fuel of an ADS. This 
result can considered as a general isation of the master equation obtained by 
Ballester and Munoz-CoboT(|2005l ) for cross statistical descriptors (problems 



with more than one detector). 

In particular, we admit the intrinsic source spontaneous fission process to 
behave as a Poissonian one, albeit neutron emission multiplicity corresponding 
to spontaneous fission events has also been included. 

In this paper we have neglected the contribution of delayed neutrons, therefore 
all quantities appearing here can be interpreted as prompt variables for short 
time-scales, in comparison with the delayed neutron precursors lifetimes. 



2 General expression for the relationship between the source pgf 
and the kernel pgf 



The derivation of the general relationship between the sour ce pgf and t he 
kernel pgf shall be based on two well known results (see, e.g., lLandol l|200ah ): 



given a random variable, Z = X\ + X 2 + . . . + X n , X iy % = 1, 2, . . . , n, being 
mutually independent discrete random variables, and Gx^s) = J2j s ^ Px i= j 
being the probability generating function associated to X^, where Px i= j is 
the probability for the occurrence of the event Xi = j, J2j Px i= j — 1> then 
the probability generating function of Z can be expressed as 

G*0O = n G *0O; (2) 

i 

given a random variable Z = X\ + X 2 + . . . + Xy, X = Xi Vz and y being 
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mutually independent discrete random variables, then 



G z (s)=G y (G x {s)). (3) 

In our system we must consider tw o indep endent sources leading to detector 
counts, Ballester and Munoz-Cobol ( 20051 ): neutrons appearing from sponta- 



neous fission of isotopes contained within the nuclear fuel, and neutrons com- 
ing from the spallation source induced by nuclear interactions of the proton 
beam with the target material. 

In particular, for the derivation of the cross correlation between two neutron 
detectors, we shall consider the two-dimensional discrete random variable Z = 
(Z\, Z 2 ), where Z iy i = 1,2, is the number of neutron detections gathered by 
the z-th detector during the time interval (t^ — r c . , tf . ) . And according to our 
previous discussion, Zi = Afi+AAi, where Afi, AAi, are the number of detections 
registered by the i-th detector and coming from the intrinsic spontaneous 
fission source and the external pulsed source, respectively; equivalently, Af = 
(Af 1 ,Af 2 ),M = (M 1 ,M 2 ). 

In case of the intrinsic spontaneous fissions source, the number of detector 
counts within the time interval (t f . — r c .,t f .), Afi, can be expressed as 

Afi = X h +X i2 + ... + X iy , (4) 



where y is the number of spontaneous fission events occurring within the 
fuel material and X^. = Xi is the number of detector counts gathered by 
the i-th detector corresponding to each spontaneo us fission event, both bein 



mutua lly independent discrete random variables, iBallester and Munoz-Cob 



( 2005h : its corresponding two-dimensional discrete random variable will be 



denoted as X = (Xi,X 2 ). We can define the following probability functions: 
P ni n 2 {d\ (tfj , d 2 (tf 2 )) is the joint probability to have Afi — Hi, i — 1, 2, detec- 
tor counts within the time interval (tf 4 — T Ci ,ifJ upon the introduction of the 
neutron intrinsic source in the remote past, P y (tf) is the probability to have 
y = y intrinsic source spontaneous fission events within the fuel material at 
time t{ = max {t^ , tf 2 } upon the introduction of this neutron source in the re- 
mote past, P XlX2 ( r , t\di (ifj , d 2 (t{ 3 )) is the joint probability to have Xi = Xi, 
i = 1,2, counts in the i-th detector within the time interval (t^ — r Ci ,tfJ 
after an intrinsic source disintegration event at time t and position r. In ad- 
dition, K nin2 (i?, t\di (tfj ,d 2 (tf 2 )) is the joint probability of having Afi = Hi, 
i = 1,2, detector counts per single neutron injected in the phase-space point 
•& = (r, v, f2) at instant t. Associated with these probabilities we have the 
source probability generating functions 

oo 

G s i(s 1 ,s 2 \d 1 (t {l ),d 2 {t {2 ))= £ s^sV?P nim {d x {t h ),d 2 (t i2 )), (5) 

71 1) "2=0 
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y=0 



(6) 



G s ,x (si, s 2 , r, t\d x (t h ) , d 2 (t&)) = ^ s* 1 s x 2 2 P XlX2 (r, t\d x (t h ) , d 2 (t h )) , (7) 



x\,X2=0 



and the kernel probability generating function 

oo 

G K (Si, S 2 , t\d ± (t f J , d 2 (t f2 )) = sT4 2 K ni n 2 t\d! (t fl ) , 4 (tfj) .( 

ni,ri2=0 



If we apply the result given by equation (3) and consider the fact that the 
spontaneous fission neut r on source will behave as a Poisso nian source, then 
we find that ([jell (|l965h : iBallester and Munoz-Cobol ^OOfih ) 



I dt J dr\ si N (t) p sf (r) [s - 1] 



(9) 



Gg (si,s 2 |di (t h ) ,d 2 (tf 2 )) = G s ,y (Gs.at (si,s 2 ,r,t|rfi (t fl ) , d 2 (tf 2 )) , tflJO) 

where we have supposed the source, which is given by the product of its time 
dependent activity, \ s fN(t), and the shape probability distribution function 
p s f (r), to be introduced in the remote past. In addition, the relationship be- 
tween the spontaneous disintegration source probability generating function 



and the kernel probability generating f unction is given by (|Munoz-Cobo et al 
( 1987b : iMunoz-Cobo and Verdul l|l987h ) 



Gs,x Oi, s 2 , r, t\di (t f J , d 2 (t h )) 
hi 



3=0 



-sf 



dv / dfi 



4ir 



G K (si, s 2 , i?, t|<ii (ifj , d 2 (tf 2 )) 



being the probability of emission of j neutrons after a spontaneous fis- 
sion within the fuel material, J s f the maximum number of spontaneous fission 
neutrons emitted after a source disintegration, and Xs{ v ) their corresponding 
spectrum. 

On the other hand, in order to obtain the relationship applicable for the 
proton-beam-driven spallation neutron source we can proceed in the follow- 
ing way: firstly, let us consider the expression for the joint probability to 



register M.i 



1,2, neutron counts in each neutron detector during 



the corresponding detector time interval (if. — T c -,tf-) following the injection 
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of one single proton (lp) belonging to one proton source pulse injected at 
a random time £ 6 (0, T) , T b eing the pulsed proton source period, i.e., 
Ballester and Munoz-Cobol ifcOOfih . 



P m \m 2 {Z\dl{th),d 2 {t h )) 



m\iri2 

mi,m2 



= J drp sp (r)E^ P E " " " E II / d ^ / d ^ 

x/sp^,^)^)^) (r^fJ^ldi^fJ,^^)), (12) 

restricted by the constraints mf + m- 2 -* + . . . m[ — i — 1, 2, where £j P is 
the probability of emission of j neutrons after a spallation interaction within 
the target material, I sp the maximum number of neutrons emitted in each 
spallation interaction, p sp (r) is the spatial distribution function for neutrons 
born after a spallation interaction within the target material, and f sp (vi, f2«) 
the spectral and angular probability distribution function corresponding to 
these spallation neutrons. If we multiply the previous expression by s^s™ 2 
and sum up from m^mj = to oo, we get the relationship between the 
spallation source probability generating function and the kernel one for one 
proton randomly injected, 



Gi p (si,S2,£K (f fl ) ,d a (tf a )) 

= / drp sp (r)X>f [T sp G K (s 1 ,s 2 ,$,Z\d 1 (t h ),d 2 (t h ))] j . (13) 

where we have used the spallation operator T sp o = J df, / df2j/ sp (v^ fij) o. 
With this expression and taking into account the relationships expressed at 
the beginning of this Section, and admitting that the number of detector 
counts after the introduction of each proton belonging to the same source 
pulse can be considered as independent discrete random variables, we can 
derive the relationship between the spallation source and the kernel probability 
generating functions for the injection of k protons in one proton pulse (pp) at 
time £, 



G pp>k (si, s 2 ,£\di (tfj , d 2 (t f2 )) 



dr Ps P (r) E £ f [Tsp G k (si, S2,&,€\di (t h ) , d 2 (tf 2 ))] J 

3=0 



(14) 



But, in general, the number of protons injected per accelerator pulse might 
be considered to be a discrete random variable, thus, the correct expression 
for the relationship between the spallation source and the kernel pgfs for the 
injection of one proton pulse can be recast as 
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Gpp (si, s 2 , CMi (*fj , d 2 (t h )) = ^ £ l PG P P ,k (si, s 2 , £\di (t h ) , d 2 (tfj) , (15) 

fc=0 

where e pp is the probability for the accelerator to inject k protons per proton 
pulse and I pp is the maximum number of protons that can be introduced in 
the system per proton pulse. 

Considering this expression, for a periodic pulsed proton source such as that 
one given by equation (1), we will have 



G s s (si, s 2 ,£\di (t fl ) ,d 2 (tf 2 )) = n G pp ( Sl > s 2,£ + rnT\di (t il ),d 2 (tf£$) 



In general, neutron detectors will register counts coming from both, sponta- 
neous fission and spalla tion, sources which can be treated as mutually indepen- 
dent random variables (Balleste r and Munoz- Cobol ifeOOfih ). thus the probabil- 



ity generating function governing both processes will be given by the product 
of equations (10) and (16), that is, 



G s (si,s 2 ,Z\d 1 (t fl ) , d 2 (t f2 )) 

= G a i ( Sl , s 2 |d a (%) , d 2 (t f J) x Gf s 2 , (t fl ) , d 2 (i f2 )) . (17) 



This last equation expresses the relationship between the source probability 
generating function and the kernel probability generating function when we 
consider the effect of the intrinsi c spontaneous fission source and t he periodic 
pulsed spallation neutron source, Ballester and Munoz-Cobol ( 2005[ ). In partic- 



ular, for the deterministic pulsed method we just need to choose the elapsed 
time £ equal to zero, whereas, in order to apply the stochastic pulsing me- 
thod we will calculate the expected value of (17), £ being uniformly sampled 
between zero and the proton pulse period, T, i.e., 



Gs (si, s 2 \di (t fl ) ,d 2 (tf 2 )) = (G s (si, s 2 ,£\di (tf 1 ),d 2 (t h ))) ( 

T 

= J — G s (si, s 2 , £\di (t fl ) , d 2 (tfa)) ■ 
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3 The Boltzmann neutron transport equation for counting pro- 
blems from the stochastic neutron transport theory 



In this Section we derive the integro-differential equation governing the kernel 
probability generating function Gk (si, s 2 , £\di (£fj , d 2 (tf a )). First of all, we 
need to obtain an expression for the probability function K ZlZ2 (i?, t\d\ (tfj , d 2 (if 2 )) 
for both neutron detectors . It can be done u s ing a proba bility balance of 
mutually exclusive events (M unoz-Cobo et al. Then we shall 

multiply the probability balance equation of i^ zlZ2 (i?, (tfj , <i 2 (*f 2 )) by the 
factor s^Sg 2 an d then sum up fro m Zj , z?. = to oo. Next , we need to apply 
the kn own Pal's methodology fseelp^ |l95^ : lBel3 (|l96^ : lMunoz-Cobo et al ' 



;o obtain the non-linear transport integro -differential equation satis - 
fied by the kernel probability generating function. iMunoz-Cobo et al\ ((1987) : 



HG K (si, s 2 , i?, £|cZi (t fl ) , d 2 (tf 2 )) 

= 5 (G K Oi, s 2 , r, ?/, n', t\di (t h ) , d 2 (tf 2 ))) 



(19) 



where we have defined the general time-dependent transport operator 

W = -f^ + n.V-E t (r ) «)^ ) (20) 

and the non-linear kernel pgf source operator 



S (o) = Et (r, v) {cf°> (r, v, £) + Sl C^ (r, u, t) + s 2 C^ (r, v, t) 
+ (cf°> (r, „, f) + Sl Cf'°) (r, t) + s 2 Cf ^ (r, «, ^ 

X J dv' J dfi'/e ( r , O') o 



where 



x 



+ E • J ( r > «. o + *icj ' ; ( r > «. o + ' j ( r > ». *) 

dv> ldn>x^±o J 

47T 



/ 8 (r,«,n|7/,n' 



/^fiKnO forr^y Dl ,F D2 , 



(21) 



represents the probability distribution function for a neutron to exit with ve- 
locity and direction within (v', v' + dv') and (ft', fl' + dft'), respectively, after 
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a scattering event (the superscript Dj applies for the i-th detector volume) with 
an incident neutron with velocity v and direction f2, whereas the spectrum of 
neutrons emitted following a fission event is given by 



X(v) for r £ VbuVba, 
X Dl (u) for r e Vbi, 



that is, as before, for the system volume no subscript is used, while to refer to 
fissions occurring within either detector volume we will add the subscript Dj. 
I is the maximum number of neutrons produced after a fission event. 

In addition, in expression (21) we must specify the C-probabilities: 



■D 2 
ic 

J>2 



if r £ V Bl ,V B2 , 
^l-^'A^y) if re V Dl , 
(l-^A(rf 2 (tf 2 ))) if rG Vb 2 , 



(22) 



is the probability to have zero detector counts following a capture event within 
the nuclear system (V S ys) o r within one of the detectors (V Di , % — 1,2) after 
a given neutron interaction at position r and time t. S t (Sj 5 *) denotes the 
neutron total macroscopic cross section for the system (i-th detector) volume, 
S c (Sj^) is the neutron capture macroscopic cross section, and rj^ accounts for 
the i-th detector capture efficiency A (o^ (tfj) — (H (t — (t { . — r c J) — H (t — t^)) 
is the time window for the i-th detector, (t) being the characteristic or Heav- 
iside function. Similarly, for one neutron count after a capture event in the 
first detector, we have 



Cf°) (v,v,t) 



ifr^l/ Dl , 
rfr^Mdi (t h ))HveV Bl . 



(23) 



The expression corresponding to the case of one neutron count following a 
neutron capture event in the second detector can be derived in an analogous 
way. 

Next, in equation (21) we must also specify the probability to have zero counts 
in both detectors after a neutron scattering event at position r and time t, v 
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being the incident neutron velocity, i.e., 



Cf 0) (r,M) 



2k 

St 



i(l-t)J*A ifrey Dl , 
jg(l-^A(d 2 (t f2 ))) ifrGK D2 , 



(24) 



where E s (E^*) denotes the neutron scattering macroscopic cross section for 
the system (i-th detector) volume and i]® 1 is the i-th detector scattering effi- 
ciency. 

In case of one detector count, for instance, in the first detector: 



C s (1 ' 0) (r,M) 



ifr£V Dl , 
fjAK fe))ifreV Dl . 



(25) 



Finally, in equation (21) the probability of occurrence of a fission event with 
emission of j > neutrons leading to zero detector counts following a neutron 
interaction at position r and time t is given by 



C 



(0,0) 



ifr^\/ Dl ,Vb 2 
^ Dl ^(l-% Dl A(rfi(t fl ))) ifreVbx, 



(26) 



where e_j (sf l ) accounts for the probability to emit j neutrons after a fission 
event within the system (i-th detector) volume, Ef (E° l ) denotes the neutron 
fission macroscopic cross section for the system (i-th detector) volume and r/ f ' 
is the i-th detector fission efficiency. Whereas, if we consider, e.g., one detector 
count registered by the first detector, 



C 



(i.o) 



r,v,t) 



ifr£V Dl , 
n^ef^Aidx (t fl ))ifrey Dl . 



(27) 



Expression (19) must fulfil the final condition G K s 2 , (tfj , <i 2 (£f 2 )) = 

1 for t > tf = max {t fl , tf 2 }, due to the causality principle, and the boundary 
condition Gk (si, s 2 , r B , v, f2, t\di (tfj , d 2 (tf 2 )) = 1 f or n • O > 0, i.e., for neu- 
trons injected outwardly at a convex boundary ([jell (|l965l ): I iMunoz-Cobo et all 
(|l987jl ). 



10 



Further, according to Bartlett's procedure. iBartlettl ((19551 1. we can derive the 
first factorial moment of the number of detector counts per single neutron 
injected in the system at the phase-space point i? and at time t as 



Zi(#,t\d{tt)) 



d_ 

dsi 



G K (s 1 ,s 2 ,^,t\d i (t { .)) 



(2* 



,S2=1 



whereas the cross second factorial moment of the number of detector counts 
per single neutron injected can be defined as 



z^{&,t\d(t { )) 



d 2 



dsids 2 



G K (si, s 2 , 0, t\di (t fl ) , d 2 (t f2 )) 



(29) 



Sl,S2 = l 



Thus, applying the operator d/dsi\ s s =1 to the expression corresponding to 

the transport integro-differential equation satisfied by Gk (si, s 2 , t\di (tfj , d 2 (if 2 )), 

equation (19), we get 



ld lAz-S+ 



(30) 



where the time-independent adjoint transport operator L + . lMunoz-Cobo et al 
fll987h . and 



x [ff (t - fa - Tc( )) 
This magnitude will be non-zero only for r G Vd 4 and £ G (tf a — r c . , £f . 



(31) 



The solution Zi of equation (30) must fulfil the boundary condition Zi = 
for n ■ 17 > 0, on a convex boundary, and the time-reversed causality condi- 
tion, i.e., it must vanish at the end of the measurement period t Ci . Indeed, 
equation (30) reveals the nature of Zi cLS £111 adjoi nt generalised Green ' s func - 
tion driven by the adjoint importance source 5p . iMunoz-Cobo et al\ ((1987). 
Consequently, for the forward transport problem we shall write 



'\_d_ 

v dt 



-L 



Si, 



(32) 



where L is the time- independent direct transport operator, Bell and Glasstonel 
fjl979h . and 



Si (#,t)=S? (^t) + Sl P (d,t) 
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= Sf (0) x Sf (t) + 5f {■&) x S x sp (t) , (33) 

that is, the total neutron source can be expressed as the sum of the intrinsic 
and the spallation neutron sources, 

Sf (0) = u 8i N oPs{ (r) (34) 
Sf(t) = X Bl N(t)/N = X a , (35) 



where we assume the initial number of nuclei corresponding to the spontaneous 
disintegration neutron source to be constant in our time-scale, 

Si P (0) = ^pp^spPsp (r) / S p (v, Q) , (36) 

00 

ST(t)= £ <5(t-(£ + mT)), (37) 

m=— 00 

with v w = J2[ w jej, w = sf , sp, pp. 

Now, the forward neutron flux satisfies the initial condition <f)(t = —00) = 0, 
and the boundary condition = for n • f2 < 0, on a convex boundary. 
A proper choice of the corresponding boundary and final conditions for the 
adjoint function makes the associated bilinear concomittance to vanish, and, 
hence, due to the commutation relation 

(zi\Si) = (S^M) , (38) 



where Dirac's notation for the inner product is used, and where both terms 
account for the average number of detector counts during its counting interval: 
at the left hand side we have the inner product of the neutron source strength, 
Si, (neutrons emitted at a given phase-space point and time) and the number 
of counts gathered at detector i per single neutron introduced at a given phase- 
space point and time, 2$; equivalently, the right hand side term expresses the 
inner product of the effective macroscopic neutr on detection cross s e ction for 
the same detector, S^ . , and the neutron flux, 0, iMunoz-Cobo et al. 



Similarly, we can make use of (29) to find an expression for the cross sec- 
ond factorial moment of the number of detector counts per single neutron 
introduced: 

(~ J^ + W=SU. (39) 
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where the importance source for the cross second factorial moment can be 
recast as 



S t (r,v) 

= C^ (r,v,t) J dv'J dtff B (r,v,n\v',rt)*i (r,i/, (%)) 

+c<°> 1 > (r,u,i) | d^y dn'/s^^niu', n')z 2 (r,v',n',t\d 2 (t i2 )) 

+ EjCf 0) / du'/ dfi'xfri/) /47TZ! (r, ft', t\d, (%)) 



+ £jCf 1) (r,v,t) / du' / dfl' X (r,v') /4nz 2 (r,v',fl',t\d 2 (t i2 )) 

3=1 

+ Ej(J"1) E ,n) (r,v,t) f dv'J d{l' X (r,v')/4n 

j=2 ri,r 2 =0,ri+r 2 =l 

xz! (r, fi', t\d, (t fl )) z 2 (r, u', SI', t\d 2 (t f J) . (40) 



The cross second factorial moment zjz^ must satisfy the same time-reversed 
and boundary conditions as Zi, but for tf = maxjt^, £f 2 }. Hence, it can be 
viewed as a generalised adjoint function, now driven by the adjoint source, 
^DiD 2 ' that is, the product of the detector cross sections and the spectral and 
angular weighted neutron importances. Again, the commutation relation leads 
to the identity 



(^|^) = (^ lD2 |0). (41) 



In order to express the adjoint problem in terms of instantaneous detector 
counting rates, we shall divide adjoint transport equations (30) and (39) by 
r c . and r cl r C2 and then calculate the limits lim Tc .^ and lim Tci)Tc2 ^ , respectively. 
As a consequence, we shall write 

Bi- L+ )*"*' (42) 

(-;|-i + )^P5=SS lDl , (43) 



where, by definition, 

W-lO-llm *** "*<**» , (44) 
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(45) 



since, taking into account equat ions (42) and (45), Zj (tf. t — tf.) can be re- 
garded as a displacement kernel, iMunoz-Cobo et all (jl987l ) , 



-=-/-q + + + + n r z i z 2 {#,t\di (t fl ) , g? 2 (tf 2 )) , v 

ziz 2 (1?, t - tf x , t - tf 2 ) = lim , (46) 

T Cl ,r C2 j.o r Cl r C2 

S D + Dl lim g^WfaMifaj) . (47 ) 

r Cl ,r 02 iO r c ,r C2 



Next, we look for a solution to the Boltzmann neutron transport equation (32) 
for the direct flux, (tf , t) , satisfying the a-modes expansion: 

<j> (0, t) = sf (tf , t) + sp (tf, f ) = E W W (*) » («) 



where we assume the eigenf unctions (£>,- (tf) to form a c omplete basis in the 
corres ponding Hilbert space (|Bell and Glasstonel ljl97flh : ICarta and D'Ang 



(1999)). These must obey the a-eigenvalue equation 



LtPi (tf) = -±<pi (tf) . (49) 
Similarly, for the adjoint flux instantaneous rate we will have 

Zl (tf, t - t h ) = E <p£ d (*) (< - <0 , (50) 

j 

and the a-eigenvalue equation 

cv ■ 

L + ^ d (*) = (*) • (51) 



Beneath this ansatz, it is obvious that we must put (^J • (tf) = (p^ 2 • (tf) = 
( tf). The adjoint and forwar d eigenf unctions satisfy the biorthogonal rela- 
tion, |Beir^nT^Glasstoni ( 1979f ). i.e., 



(~^n> Vm) = <*nm (^n , Vr^J , (52) 



where the phase-space inner product is defined by (a, b) = J drdfdfia (tf) b (tf). 
If we introduce the ansatz (48) in equation (32) and apply the Fourier trans- 
form operator to both sides of it, we shall obtain the Fourier transform of 
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the j-th flux instataneous rate time-dependent term, which, on account of 
identities (49), (51), and (52), reads as 



iuo — a 



3 



with 



(Sf, <pf) = I dr J dv J dnu s{ N oPsi (r) (r, «, O) , 

(Si P , V 3 /) = J dr J dv J driz/ pp z/ sp p sp (r) / sp (v, n) <£>+ (r, v, fi) 



(54) 



(55) 



In the same way, from (42) we can deduce the expression corresponding to the 
Fourier transform of the j-th adjoint flux instantaneous rate time-dependent 
term for the z-th neutron detector: 



VjJ + a 



(56) 



where S£. = S£. (•&) is given by (31). 



4 Analytical expression for the cross power spectral density with 
pulsed sources 



4-1 The deterministic pulsing method 



In order to obtain the analytical expressions corresponding to the factorial 
moments of the number of counts of both neutron detectors we can apply 
again Bartlett's procedure ( Bartlett ( 1955h ) to the expression corresponding 
to the source probability generating function. As we have outlined previously, 
we can do it taking into account two different situations: in the first case, 
we can calculate factorial moments corresponding to the deterministic pulsing 
method. 



We are interested in the well known cross covariance function, Papoulisl (1991) 
defined as 



(di (t fl ) ,d 2 (tf 2 )) 



d 2 G, 



dsids 2 



Sl,S2 = l 



dsi 



dG, 



x 



Sl,S 2 =l 



8sq 



(57) 



Sl,S2 = l 



15 



where Gs = Cs (si, s 2 \di (%) , ^2 (*f 2 )) is given by (17) with £ = 0, i.e., it 
is the difference between the cross second factorial moment of the number 
of detector counts gathered by both detectors and the product of their first 
factorial moments. It can be recast as 

S (di (t h ) , d 2 (t f J) = (S+ D J0) (dx (t fl ) , d 2 (f fa )) + AS (d 1 (%) , d 2 (£ f2 )) , (58) 



where the first term takes into account the contribution coming from multi- 
plicative processes within the system and the detector volumes due to fission 
events and detections. The latter can be de facto neglected if we admit that the 
volume occupied by detectors is small in comparison with the system volume. 
In addition, the second term in (58) stems from the non-Poissonian behaviour 
of both neutron sources. 

Next we can divide equation (58) by r Cl r C2 and then apply the limits lim Tci Tc2 j_ 
in order to derive the expression corresponding to the second order instanta- 
neous rate of the cross covariance function: 

E(t-t h ,t-t h ) = (S+ lB2 \<f>)(t-t h ,t-t h ) + AZ(t-t h ,t-t h ), (59) 

where, by definition, 

E(t-t„ t - t[2 )= lim 5(*fa>.*fa» , (60) 

= ito m<hM,d,M) i {6l) 

T Cl ,Tc 2 J.U T Cl T C2 



with 



(S^dM (t-t h ,t-t h ) 

= (^ lD j0 sf ) (* - t h ,t- t h ) + (£+ D2 |<F> (t - t h ,t- t i2 ) 

tf 

D = v{y — V)jv 2 being system Diven's factor, v = Y,[j £ j, v {y — ^) = 
X) 2 3 {j ~ an d where we have defined the phase-space inner product 



(pt<Pj, <Pk, <pf) = j dr j dv j dfi£ f (r, v) ifj (r, v, Q) 

Vsys+Vd 
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x / dv" I dn"^^-tf (r, v", Q") , (63) 

and 

AS (t-t h ,t- t h ) = A5 sf (t - t h , t - t fe ) + AH sp (t-t h ,t- f fa ) , (64) 

which arises from the non-Poissonian nature of the intrinsic and the external 
spallation sources, respectively, 

AS sf (t-t h ,t-t i2 ) 

t f 

= £ (Sf, tf, tf) J dt&j (t - t h ) (+ 2k (t -t i2 ), (65) 



AH sp (t-t h ,t-t i2 ) 

= e M» (s?, tf, tf) + p PP - 1) Or, tf) Or, vt)_ 

3,k 
tf oo 

x f dt £ Ht-mT)C£ l3 (t-t h )(+ 2k (t-t i2 ), (66) 

-oo m=-oo 

where D w = z/ w (z/ w — \)jv\, w = sf, sp,pp, is Diven's factor for the sponta- 
neous fission (intrinsic) source, the spallation neutron production source, and 
the pulsed proton source, respectively, with z/ w (z/ w — 1) = ~~ ^) £ J- Fur- 

thermore, in the last two expressions we have introduced the following inner 
products: 

(sf, tf, tf) = u a N J drp sf (r) 

x Jdv'jdn'^l^(r,v',n') 

x jdv" jdn"*^tf{v,v",Sl"), (67) 



(ST, tf, tf) = VppVsp j drp sp (r) 

xjdv'j dn'f sp (v',n')tf(r,v',n') 

x / dv" j dQ"f sp (v", Q") tf (r, v", Q") . (68) 
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Without loss of generality, we can assume that the upper integral limit in (57) 
can be selected in such a way that tf > max {t^, tf 2 }, and then admit that 
tf — > oo. Next, we can define the time delay between the final instant of both 
detector intervals as r = t f2 — t fl , and then apply the operator / drexp (— iur) 
to equation (59) to derive the Fourier transform of S (t — tf i: t — tf 2 ), i.e., the 
cross power spectral density: 



CPSD = ^ (S£ lD » +F 



AS 



(69) 



where, on account of equations (53) and (56) for £ = 0, 



^ 1D2 I0 SP 



(70) 



S 



DiD 2 l 



isf 



X 



A 



sf 



(71) 



S + I<M1 - V u 2 D fee ^ ^) (^ + ) (gjjjgO (g^jgO 

3,k,l [v ( Pj>'Pl){v'Pk,<Pk)[v<P7,<Pl) 



X 



[(exp(-a/T)-l) 1 -(exp(-(a fc + i^)T)-l) 1 



(a; + i (a, - a k )) (u + ia z ) 



where, owing to the application of the deterministic pulsing method, we just 
need to assume that t tl = I P T, with I P — > 00, whereas 



AS 



AS sf 



AS sp 



(73) 



A^sfl _ V -2 r, ( csf ( .+ ^ (ggiljgj) ( g D2^fc 
3,k {vVj'Vj) {vVkiVk) 



x 



[uj - \ai 3 ) {uj + ia k y 



(74) 



j,k 
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X 



(■Spi'lj) ( g D 2 ^fc) (exp (- (aj + ico) T) - 1) 



-i 



(75) 



4-2 The stochastic pulsing method 



Similarly, expressions for the factorial moments of the detector number of 
counts can be derived for the stochastic pulsing method, just applying Bartlett's 
procedure to expression (18). Now the cross covariance function will be defined 

as 



(E(d t (t h ),d 2 (t f2 ))>. 



d 2 G, 



dsids 2 
/dG s 



Sl,S 2 = l / 



ds i 



/ dG, 



Sl,S 2 =l/ £ 



\ ds. 



(76) 



,S2=1' 



where Gs = Gs (si, S2, (tfj , c?2 (*f a )) * s given by (17), which can be recast 
as 



(S (dx (t fl ) , d 2 (t f2 ))) f = «%|^ p > <^ p » 5 - «^ p » 5 

+ (z^\S 1 ) ^ + (AE) V (77) 

where S^ 13 = Sl P t) is given by the second term of expression (33). The 

first and second term of (77) s tem from the time correl a tion introduced by the 

stocha stic pulsing method (see lCeder and Pazsit ( 2003| ) : Ballester and Munoz-Cobol 
( 2005h ) and, obviously it only involves the pulsed neutron source. 



The cross power spectral density shall be derived following the same steps as 
before, i.e., dividing equation (77) by t Ci t C2 , then applying the limits lim Tci jTc2 1 
to the expression obtained, and, finally, calculating its Fourier transform. 
Thus, we shall write, on account of the commutation relations (38), (41), 



(CPSD) £ = J^ 



DiD 2 l 



AS 



Di 



D 3 I 



(7S 



where, as before, 



T 



'(SUM) 



T 



DiD 2 i 



(79) 
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T 





= T 


AS sf " 















(80) 



The external pulsed source contribution to the term arising from the system 
cross covariance, equation (79), under the stochastic pulsing method becomes 



T 



D 2 I 



X 



J2" 2 D \pi<Pj><Pt> ( Ptj JT 



Vk ) ^) (-«,) (w - it**) (w + ia,) 



On the other hand, the time correlation term stemming from the stochastic 
pulsing method in (7 7), after dividing by r^r^ and c alculating the limits 



lim^o, is gwen bv. iRallester and Miifinz-Cnhr] J^OOsh. 



j,k 



-oo — oo 



x&^-^-r) E expfi— (f-f; 

m=— oo,m^:0 



2irm 



(82) 



and applying now the operator / drexp (— iur), we find that 



2?r (Sp^Vj) (gSa^fc 



CJ - —J 



/„ (w- ia,) (a; + iaa) 



We can similarly proceed to calculate the second term appearing in equation 
(80): 



E 



[u sp D sp (ST, tf } <p+) + (D pp - 1) (S?, (S?, ^ 
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X 



1 



T (uj — ia 7 ) (a; + ia^) 



J4) 



This term can effectiv ely become negative for certain experimental conditions, 
Ballester et~aJ\ (|2005[) . 



5 Fundamental mode approximation and discussion 



Following the derivation obtained in the previous Section, we can also apply 
the fundamental mode approach to the expression corresponding to the CPSD 
when the stochastic pulsing method is used: 

( cpsD)^( ft+fe t X"-T))whs>- (85) 



Ox = I v 2 D (Sfv? , <Po, <Po) -77 



sf+sp 



¥>d 1 



(-«o) 



+ 



1 r 

T 



Q2 



(uVo^o) (J^o^o)' 

/V p t»+\ 2n ( s ^ ,<Po ) ( s ^ Vo ) 



(86) 
(87) 



where we have defined the total neutron source strength S{ +sp = +sp (i?, t) 



Equat ion (85) apparently seems to be similar to that obtained bv lRugama et al. 
( 2004h albeit a deeper examination permits to understand an important dif- 
ference: in the latter case, the expression for the CPSD for the stochastic puls- 
ing method was derived assuming a quasi-Poissonian behaviour of the pulsed 
external source, i.e., with no delay time averaging. That means that the ex- 
pression obtained in this reference should be exactly equal to our expression 
for the CPSD in Section 4 under the deterministic pulsing method, only with- 
out the contribution of (75), which stems from the non-Poissonian behaviour 
of the periodic pulsed source. But it is clear that, in any case, equation (69) 
is a completely bounded function for any non-zero value of the frequency, u, 
which cannot produce the response obtained at MUSE-4 experiments, also 
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-ato) obtained for the configuration SCO (with pilot rod inserted) of 



Table 1 
Values of (- 

the M ASURCA subcritical assembly used during MUSE-4 experiments. ISoule et al. 
(2Q0J). The result of the eigenval ue fitting by means of the Feynman-a method is 



reported bv lBallester et al 
Method 



( 20051 ). 



-Qo) (rad 



Stochastic Pulsing CPSD method 13258±273 
Stochastic Pulsing Feynman-a method 13646±515 



reported by iR.ugama et al\ ([20041 ) . for the stochastic pulsing CPSD method. 
Unlikely, we have shown that those spectral lines appearing at frequencies 
which are multiples of the accelerator frequency are indeed not produced by 
the system cross covariance contribution, but they are a prima facie of the 
time self-correlation introduced by the stochastic pulsing method. 



In addition, it has been recently shown that the utilisation of determinis- 
tic external pulsed sources, such as those used in MUSE-4 e xperiments, can 



make the nuclear process to behave as a sub-Poissonian one. iBallester et al 
( 200,51 ). [t can occur when the contribution of non-Poissonian term of the 



pulsed source becomes negative. That means that under those conditions the 
time-independent term Qi can be effectively negative for very deterministic 
pulsed sources. But, anyway, it is of common practice to use Bode's diagram 
of the magnitude in order to represent graphically a system response function 
such as (85), therefore, it is not necessary to consider the sign of Qi provided 
that, from a practical viewpoint, we only need to fit its magnitude for those 
points of the graph where uo is not a multiple of the accelerator frequency, 
together with the eigenvalue «o- Notice that the value of Q2 cannot be de- 
termined from this technique (we will need a further integral condition), but 
in any case it is completely useless for a practical purpose. That means that 
when the stochastic pulsing CPSD method is applied to determine the value 
of the subcriticality level of a nuclear system, only two fitting parameters must 
be considered, in contrast with the sto chastic pulsing Feynm an-a; technique, 



where three parameters must be fitted, IBallester et al. (2005). 



In Figure 1 we show some experimental points report ed in that graph cor- 
responding to the stochastic pulsing CPSD method of iRugama et al\ (|2004f) 
(Figure 6), together with equation (85) conveniently fitted. This particular 
experiment, corresponds to the configuration SCO of the MASU RCA subcrit- 

" (|2004t ). with a D-D 



Soule et al 



ical assembly used during MUSE-4 studies, 
pulsed source. The v alue of the prompt neutron time constant obtained by 
Rueama et al\ ([20041 ) . which is shown in Table 1, is effectively very similar to 



that reported using other noise techniques for the same conditions. 
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Fig. 1. Stochastic CPSD pulsing method obtained during MUSE-4 experiments 
for the SCO configuration of th e MASURCA subcrit ical assembly. Points of the 
experimental curve (points) from R ueama et al\ ( 20041 ') are compared with equation 
(85) fitted with (— ao) = 13258 ± 273 rad • s _1 (solid). The accelerator period is 
equal to 1 ms. 



6 Conclusions 



In the present work we have dealt with the applicability of stochastic- neutron- 
field-based methods for the study of the neutron counting statistics in a nuclear 
system. We have derived the generalised two-detectors relationship between 
the probability generating functions of the kernel and the source for sub- 
critical assemblies when pulsed neutron sources are used together with the 
intrinsic neut ron source coming from spontane ous fission events within the 
fuel material, Ballester and Munoz-Cobol ( 20051 ) . It has been done within the 



stochastic neutron transport theory framework, which permits to understand 
how the general transport problem is influenced by its spatial, spectral and 
angular dependence. 

Further, we have followed Pal-Bell's methodology for the derivation of the 
integro-differ ential Boltzmann tr a nspor t equation, and applied the formalism 



described by iMunoz-Cobo et al\ ()1987t ) in order to calculate the chosen sta- 



tistical descriptor. 

In Section 4 an expansion in a-eigenvalues for the cross covariance and the 
CPSD of two-detectors stochastic counting rates have been obtained. The 
contribution of higher harmonics in subcritical monitoring problems shall play 
an important role in ADS assemblies. In this case, the excitement of higher 
modes could be relevant in situations of normal operation, and it will increase 
as the reactor departs from the criticality condition. 

In Section 5 we have compared the expression obtained for the stochastic 
pulsing CPSD method with experimental data obtained during the MUSE-4 
European project. The value of the prompt neutron time constant fitted is 
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comparable with others methods. In addition, the reduced number of fitting 
parameters makes this method suitable as a subcriticality monitoring tech- 
nique. 
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